LL 


security  classification  or  this  page  fwht*  o*,  £»•»« 


hej:ffects  of^ross-^lcJTIUWal  \ 

JjE  LOCATION  ON  JRAPPED ELECTRON/ 


/ 


) REPORT  DOCUMENTATION  f^AGE 

’■■■“ — “ K/'1-  GOVT  ACCESSION  NO; 


‘Qj/RL'MK-3L±r 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


*.  PERFORMING  ORO 


~ ""  I »•  COmVraCT  OR  GRANT  NUMAERf^ 

OROAN?IATI0^3e  ANO  vL^1  ~I0.  PROORAM  ELEMENT.  PROJECT.  TASK 


I.  PERFORMING  ORGANIZATION  NAAIE 

Naval  Research  Laboratory 
Washington,  D.  C.  20375 


AREA  A WORK  UNIT  NUMEERS  ' 

NRL  Problem  H02-37 
Project  ERDA  E(49-20j 
1006) 


<1.  CONTROLLING  OFFICE  NAME  ANO  AOORESS 

Energy  Research  and  Development  Admin 
Washington,  D.  C.  20545  I 

IT  monitoring  agenCv  name  a aoor2sJ777 SmSSm  S5S  Caumiiti  ottun  | u.  security 


UNCLASSIFIED 

oec^?f,cat,on/oo-ngrao,hg 


Approved  for  public  release;  distribution  unlimited. 


17.  OtSTfllSuTiON  STATEMENT  (’•I  tfia  abstract  an  f rad  In  Block  20,  It  dlltarant  Irom  I f abort) 


t#  supplementary  notes 


This  research  was  sponsored  by  Energy  Research  and  Development 
Administration  under  Subtask  (49-20:1006). 


19.  KEY  WOflOS  (Continue  on  raaataa  atda  It  naoaaaory  and  tdaotlty  by  block  nttmbar) 


Tokamaks 

Trapped  electron  modes 
Cross-section 
longation 


20.  «ST9ACT  (Conilmto  on  raaaraa  a! da  It  nacaaaory  and  tdaotlty  by  block  nvmbar) 

"^The  dispersion  relation  for  the  trapped  electron  mode  in  non- 
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THE  EFFECTS  OF  CROSS-SECTIONAL  ELONGATION 
ON  TRAPPED  ELECTRON  MODES 


I.  INTRODUCTION 


Recently  It  has  been  shown1  that  the  toroidal  drift  of 
trapped  electrons  can  have  a large  destabilizing  effect  on 
the  dissipative  trapped  electron  mode.  This  Is  particularly 
true  for  high  temperature  regimes  where  the  dissipation  due 
to  scattering  of  the  trapped  electrons  becomes  small.  The 
destabilization  may  be  viewed  as  resulting  from  a wave-particle 
resonance  occur ing  when  the  toroidal  component  of  the  wave 
phase  velocity  is  equal  to  the  trapped  particle  toroidal 
drift.  On  the  other  hand,  Glasser  et  al.2  have  shown  that  the 


cross-sectional  elongation  of  a tokamak  can  slow  down  or  re- 
verse the  toroidal  trapped  particle  drift.  Here  we  formulate 
and  evaluate  the  relevant  dispersion  relation  to  study  these 


Note:  Menuecript  submitted  September  IS,  1977, 


II.  COORDINATES 


To  begin,  it  is  necessary  to  define  the  coordinate  system 
Let  ¥ denote  the  poloidal  magnetic  flux  function,  such  that 
surfaces  of  constant  t|)  are  magnetic  flux  surfaces,  and 
dty  - RBgdl^,  where  BQ  is  the  poloidal  component  of  B,  R is  the 
tokamak  major  radius  coordinate,  and  dl^  is  a differential 
length  perpendicular  to  the  flux  surface.  Also  let  6 be  a 
coordinate  denoting  the  poloidal  position  on  a flux  surface. 


q(  )RBa  ( 9 ) 


toroidal  component  of  B,  dlfl  is  a differential  length  along 


the  line  of  constant  in  the  poloidal  plane,  q(ty)  ■ <211 ) 

J*  ldl0'  and  the  subscript  c on  the  integral  symbol 

c #. 

denotes  integration  over  one  circuit  around  the  flux  surface 


Note  that  RB,.  is  constant  on  a flux  surface.  Ve  assume  that 


the  mode  is  localized  to  a mode  rational  surface  q(ty)  5 m/n 


III.  DISPERSION  RELATION 


Ve  now  turn  to  the  formulation  of  the  dispersion  re 


lation.  Ve  assume  that  the  wave  is  purely  electrostatic  and 
neglect  effects  nonlocal  in  (a  treatment  of  shear  stabiliza- 
tion in  non-circular  cross-section  appears  in  Ref.  3).  The 
perturbed  distribution  function  Fj  for  species  j is  given  by 
the  linearized  Vlasov  equation, 


where  $ is  the  perturbed  electrostatic  potential,  B is  the 

m0 

equilibrium  magnetic  field,  and  Fj,  the  equilibrium  distribution 
function,  is  given  by 

Ts  - *5  [x  - fBa  - • <s> 

V ' ' 'r  *'  <4-  ’ < i j * 

In  Eq.  (2),  is  the  Maxwellian  distribution  function, 

Vj  = (2Tj/nij)i  , Va'X  S - d/n  n0/d¥,  nj  5 [d£n  Tj/d*]/ [din  nQ/df]  , 

and  Tj,  nQ  are,  respectively,  the  temperature  and  unperturbed 
electron  or  ion  density. 

Since  we  expect  the  wave  electric  field  to  be  nearly  per- 
pendicular to  B,  we  express  the  perturb at; 3n  in  the  form 


£ « $(9)  exp(im9  - in  / - ieet)  (3] 

where  we  expect  $(0)  to  be  a s'i-.'\  function  of  0 compared  with 
the  variation  in  the  exponential  part  of  Eq.  (3).  Since  we 
employ  the  local  approximation,  any  Y dependence  of  <J>  is  • 
neglected  by  Eq.  (3). 


a.  Ion  Response 

To  evaluate  the  ion  density  response,  we  substitute  Eqs.  (2) 
and  (3)  into  Eq.  (1)  and  Integrate  along  unperturbed  particle 
orbits, 


im0(t'  ) - in;(t' ) - iut 


we  write  the  ion  orbit  as 


(t )sin6 


where  6 is  the  angle  of  v with  respect  to  the  direction  perpen 

dicular  to  the  flux  surface,  is  the  ion  cyclotron  frequency 

2 2 

and  we  have  assumed  that  >>B(,  and  a large  aspect  ratio  so 
that  ft.  and  B are  approximately  constant  on  a flux  surface. 


Also  we  have  assumed  the  scale  length  for  changes  in  Br  perpen 


dicular  to  a flux  surface  is  long  compared  to  an  ion  Larmor 


radius.  Putting  these  orbits  into  Eq.  (4)  and  making  use  of 

n 

the  identity  exp(iY  sin  a)  Jn(Y)^xp(  inoi) , 

ns— co 

whore  J.,  is  Bessel  function  of  order  n,  we  obtain  for  w«Q. 


tme-in? 


-iwt  * 


(8) 


Note  that  the  argument  of  the  first  Bessel  function  and  v are  to 
be  evaluated  following  the  orbit  while  the  argument  of  the  second 
Bessel  function  is  evaluated  at  the  Eulerian  position.  Since 


A 

both  $ and  JQ  a*e  not  rapid  functions  of  0 we  may  neglect  0^  in 


0(t’)  appearing  in  Eq.  (8).  To  evaluate  the  ion  density  we 
integrate  Eq.  (8)  over  all  velocities.  Interchanging  the  order 
of  the  t’  and  integrations  yields 


5n 


i e<P 

“ — - sr- 


»r 


,e 


im0  + 


i^-e 


inC  f 


expC-v^  / Vi2)  f t 


dv 

J ° 

J —00 


T 

ir*v4 


I! 


~iut  * 


dt'e 


$ C01) 


‘ w. 

U)  + 1 


i-Oj  + n4 

3 1 1 


)]  H(o,a'  ) + M*\  G(ata') 

' J 2T 


where 


H(a,a' ) = exp  |-(a2  + a'2)/4j  I0(ao'/2)  , 

G(a , o'  ) = exp  [-(a2  + a' 2)/4|  [act' 1^(00' /2) 
- (a2  + a'2)  I0(aaV/2)/2]  , 


5 


In  obtaining  the  above,  we  have  done  the  v 
use  of  the  integral 


xdx  exp(**p  x )J0(ox)J0(8x)  =(2p2)  exp[.(o2+  82)/4p2 
I <a8/2p2). 


a exp(ip0 ' ) 


The  t'  integration  is  now  easily  done.  Finally  the  v integration 
may  be  expressed  in  terms  of  the  plasma  dispersion  function  Z(5)> 


sjr  exp(im0-in5-iut)£  *p.  «pp 

1 P.P' 


exp(ip0) 


*Z(C> 


(12) 


where  w*e  5 w*  (l-|ne  +nev2/ve2),  veff  2 E-1ve(v8/v)3 

is  the  effective  collision  frequency  for  trapped  electrons,  e is 

the  inverse  aspect  ratio  (e  5 _) ,va  is  the  electron 

mm'  max  e 

collision  deflection  frequency,  denotes  the  velocity  space 
integral  over  the  trapped  electrons,  angle  brackets  denote  bounce 
average,  i.e.  <$(0)>  - Tb*"1fdle(B/B6)v|1*1  3(e),  tb  £ fdlQ(B/B0)  v^1 

i 12  —1 

v « v(l-XB)  , X £ uC^m^v  ) specifies  the  velocity  pitch  angle 
of  an  electron,  u is  the  magnetic  moment,  § denotes  integration 
over  the  electron  bounce  orbit  whose  limits  are  defined  by  XB  » 1 
J = §Vf\  (B/Bq)  dl(  is  the  longitudinal  adiabatic  invariant,  and 

ud(X,v)  £ nmec(etb)"13J/3i{(.  (13) 

Again  using  the  expansion  in  Eq.(10),  we  obtain 


6n 


e _ e 


exp(im0-in;-iut)]£)  a ,N  , exp(ip0), 

x**  PfP’  v 


where 


2ir  ■»„{» 

- /*d0e  f d3’ 


(“-“*e)FS  <exPCiP,0)> 


1 


(14) 


c.  Dispersion  Relation 


Applying  the  quasi-neutrality  condition  we  obtain  the  matrix 
equation 

« r ^ 

(15) 


z 

p* 


™pp-<"> + Npp-  <“> 


ap,  - 0. 


The  dispersion  relation  is 

det  [TMpp,(u0  + Npp'  (U>)]  * °*  ’ (16) 

Once  m is  obtained  from  (16)  thepololdal  mode  structure  can  be 
found  from  (15).  Of  course,  in  practice,  one  must  do  this  by 
truncating  the  infinite  matrix  to  a finite  matrix.  This  program 
has  been  carried  out  for  a specific  model  equilibrium  (described 
in  Sec.  IV)  and  the  results  are  reported  in  Sec.  V. 


IV.  EQUILIBRIUM 

To  evaluate  the  various  quantities  involved  in  the  dis- 
persion relation  and,  in  particular,  the  toroidal  drift  frequency 
(Dg,  we  first  establish  a solution  to  the  magnetohydrodynamic 
equilibrium  equation: 

(fp  - 5 k * ^s)*  - - 4tb2  I?  » - H*  <»V2'  <17> 

where  z is  the  vertical  coordinate  normal  to  the  toroidal  plane 
and  p is  the  plasma  pressure.  Following  Reference  2,  we  specialize 
to  a particular  analytical  model  equilibrium,  namely,  an 
equilibrium  in  whch  the  magnetic  surfaces  are  nested  ellipses 
of  the  same  ellipticity,  < , and  the  to-oidpl  current  density 
is  constant  (i.e.,  the  right  hand  side  of  Eq.  (17)  is  constant. 
Letting  x ■ R-RQ  « R0.  Eq.  (17)  reduces  to 

/a2  a2  \ 

\ — ^ i|>  * constant.  (18) 

*3x  3z  / 

Equation  (18)  yields  the  solution 

«*  - <BQp2/2q,  (19) 

where  x * p cos  6 and  z ■ <p  sin©,  B(  is  the  magnetic  field  at 
p * 0 and  q is  constant.  It  is  easy  to  show  that  the  coordinate  6 


10 


defined  in  Sec.  Ill  is  the  same  as  the  elliptical  parameter  9 
defined  above. 


The  bounce  period  tb  and  the  longitudinal  adiabatic 
invariant  J can  now  be  written: 

Tb  * KB0v)"1f d9  (20) 

J - qvBQ  1j  d0  RB(1-X‘B)*,  (21) 

2 2 

- e cos9  + e cos  0 
, and  e 2 p/RQ. 

The  reflection  positions  (i.e.,  the  points  where  v#  - 0),  Bc, 
for  a trapped  particle  with  pitch  angle  variable  X is  determined 
by  the  equation 

XB(0c)  - 1.  (22) 

If  ic  is  sufficiently  large  such  that  2e  + (<2-l)  eq”2>l, 
the  trapped  particles  will  assume  two  types  of  orbit: 

(i)  If  XB(9  « 0)<1,  Eq.  (22)  allows  only  one  positive 
solution  for  9 (<tt),  hence  the  trapped  particles  bounce  between 

C * ....... 

-0  and  0 , as  in  the  case  of  circular  cross  section  tokamaks. 

c c 


where  R»RQ  + p cos0,  BsBQ  j ] 
+ | e2q"2[l  + (k2-1)  cos26]| 


(ii)  If  AB(0  * 0)>1,  Eq.  (19)  allows  two  positive 
solutions,  0^*^  and  » and  the  trapped  particles  will  bounce 
between  and  0^2\  or  between  -0^^  and  -0c**^* 

Substituting  Eqs.  (19)~(22)  into  Eq.(13)  and  expanding 
in  terms  of  c,  we  obtain 


where 


nm 

“d-  — e 2 h (X,q,c) 
eB0p2 


h(Afq,e)  S - /*6  *i/X. 

3 


xx  S 1 XBoqic"1|-eq“1cos0  + 2e2cos20+  e2<T2  [n-(ic2-l)coe20]J 

“ c2tc”V1(l-AB)  [l  + (k2-1)cos2o]  j l«>|e2q^l^j(KS-l)eos8e]j  , 
X2  s (1-XB)*  j l-|e2q"2  [l  +( k2-1)c6s2o] J, 


X3  s Jl  + 1 e2q”2  [l  + (ic2-l)cos2e 


\ 

l]  | (1-AB)'*. 


Figure  1 shows  h as  a function  of  XBmin  for  q ■ 2, 
e * 0.25,  and  several  values  o'f  <,  where  Bm^n  is  the  minimum 
magnetic  field  on  the  flux  surface  e « 0.25.  Note  that  for 


ic>  3.5,  all  trapped  electrons  have  negative  toroidal  drift. 


V.  RESULTS 


As  shown  in  Eq.  (16),  the  full  dispersion  relation  is  the 
determinant  of  an  infinite  matrix,  each  element  of  which  is 
an  infinite  sum  resulting  from  the  ion  orbits.  This  together 
with  the  fact  that  the  growth  rate  y is  generally  comparable  to 
the  frequency  «r  renders  the  usual  perturbation  techniques 
useless  for  obtaining  an  approximate  analytical  solution  of 
Eq.  (16).  We  have  therefore  numerically  evaluated  Eq.  (16) 
for  the  model  equilibrium  of  Sec.  IV.  The  infinite  sum  in 
[Eq.  (12)]  is  evaluated  by  including  aa  many  terms  as  needed. 
Excellent  convergence  is  always  obtained  by  truncating  the  sum 
at  |s(<10.  Finally  the  determinant  is  truncated  at  the  point 
of  converge dee  (typically  a 7x7  matrix  is  sufficient)  and  evaluated. 
Two  checks  have  been  made  to  insure  that  the  mode  with  maximum 
growth  rate  has  not  been  left  out.  First,  the  Nyquist  technique 
was  used  to  monitor  the  number  of  roots.  Secondly,  in  the  limit 
of  low  mode  number  £a<<l,  see  Eq.  (9)J  , the  matrix  in  Eq.  (16) 
is  essentially  diagonal  and  the  root  reduces  to  the  root  of 
the  first  diagonal  element  as  expected.  Since  y<<wr  in  this 
limit,  the  numerical  value  of  the  root  can  also* be  checked  against 
approximate  analytical  expressions. 


In  what  follows,  we  assume  a hydrogen  plasma  with  Te  ■ T^, 
» -'i  » 1,  s = 0/R  - 0.25,  q ■ 2,  and  Z ff  - 2 . The  density 


scale  length  L = -n  (dn./dp)”1,  is  specified  as  L. 


0.2R„ 


The  parameters  BQ , nQ , and  Te  can  be  scaled  out  of  Eq.  (16),  if 


13 


we  normalize  w , , and  vef f with  respect  to  the  electron  bounce 

frequency  (“>be)  and  treat  v*  (s  veff/u>be)  and  nqp^p  as 

variable  parameters,  where  - vimic/eB()  and  n is  the  toroidal 

mode  number.  However,  to  be  specific,  we  let  BQ  ■ 45  kG, 

13  -3 

ne  ■ 5x10  cm  and  treat  T#  and  n as  variable  parameters 
instead.  {The  data  so  obtained  can  be  converted  into 
the  more  general  representation  as  described  aoove.)  With  these 
parameters  specified,  we  now  solve  Eq.  (16)  for  * (■  «r  + iy), 
varying  the  parameters  *,  Tft,  and  b.  Figure  2 shows  plots  of 
the  growth  rate  maximized  over  mode  number  versus  electron 
temperature  for  several  different  values  of  *.  Figures  3a  and  b, 
which  display  »r  and  y as  a function  of  n for  T#  • « 3 keV, 

are  typical  dispersion  curves.  From  Fig.  2,  we  note  that  for 
* • 2,  3,  and  4,  the  maximum  growth  rates  are  reduced  by  factors 
of  0.8,  0.6,  and  0.4,  respectively,  as  compared  to  the  circular 
case  (*  * 1),  and  this  is  approximately  independent  of  electron 
temperature.  Although  these  growth  reductions  are  modest, 
they  may  still  be  significant  since  they  imply  that  the  amount 
of  shear  necessary  to  stabilize  the  mode  is  correspondingly  re- 

Q 

duced.  Furthermore,  Fig.  3b  shows  that  the  mode  number  for 

which  y peaks  increases  as  the  ellipticity  increases.  This 

might  be  significant  if,  as  is  sometimes  argued  (e.g.  turbulent 
2 

diffusion  'V'  y/k  ) , the  short  wavelength  modes  are  less  dangerous 
than  the  long  wavelength  modes.  We  note  that  although  the 
electron  toroidal  drift  has  been  reversed  for  all  electrons 


for  tc  ■ 4,  there  is  still  no  dramatic  stabilization  of  the 
mode.  This  might  be  because  for  tc  * 4 there  exist  modes 
with  reversed  toridal  phase  velocity  (cf.  Fig.  3 which  shows 
a region  of  negative  wr  for  n > 290.) 

7 

It  is  interesting  to  note  that  some  equilibrium  studies 
8 how  that  even  if  the  elongation  of  the  plasma  boundary  is 
modest,  the  interior  magnetic  surfaces  can  become  very  elongated. 

One  of  us  (KRC)  would  like  to  thank  Dr.  A.  T.  Drobot, 

Dr.  I Haber,  Dr.  W.  E.  Hobbs,  and  Dr.  V.  Jones  for  helpful 
discussions  on  the  numerical  aspect  of  the  problem.  This  work 
was  supported  by  the  U.S.  Energy  Research  and  Development 
Administration . 


Fig.  3 - wr  and  y versus  toroidal  mode  number  n for  T,  - Tj  • 3 keV. 
Other  parameters  are  as  indicated  in  the  caption  of  Fig.  2. 


1 1) 

'• 

¥ \ 

3 

___ 

To  2 

:! 

S 

n 

¥ ‘ l 

Q 

Ww  1 

■ 

3 

REFERENCES 


1.  J.  C.  Adam,  W.  M.  Tans  and  P.  H.  Rutherford,  Phys.  Fluids 
19.  561  (1976),  also  W.  M.  Tang,  P.  H.  Rutherford, 

H.  P.  Furth,  and  J.  C.  Adam,  Phys.  Rev.  Lett.  35,  660  (1975). 

2.  A.  H.  Glasser,  E.  A.  Freiman  and  S,  Yoshikawa,  Phys.  Fluids 
17,  181  (1974). 

3.  E.  Ott,  W.  M.  Manheiraer,  and  K,  R.  Chu  (to  be  published). 

4.  B.  B.  Kadomtsev  and  0.  P.  Pogutse,  Doklady  Akad.  Nauk  SSSR 
186.  553  (1969)  [Soviet  Phys.  - Doklady  14.  470  (1969)]. 

5.  L.  Spitzer,  Jr.,  Physics  of  Fully  Ionized  Oases  (Interscience, 

New  York,  1962)  p.  132. 

6.  It  is  shown  in  Ref.  3 that  shear  induced  damping  rate  is 
almost  independent  of  ellipticity. 

7.  J,  P.  Friedberg,  Bull.  Am.  Phys.  Soc.  21,  1033  (1976). 

■ m 

I 1 
| 


TOROIDAL  CTR  DISTRIBUTION  LIST 


UNCLASSIFIED 


Oak  Ridge  National  Laboratory 
P.  0.  Box  X 

Oak  Ridge.  Tennessee  37830 

Aten:  J.  Clarice 
J.  Callen 
J.  Hogan 
D.  McAlees 
M.  Roberts 

la  Ra  Da  A* 

Germantown,  Maryland  20767 

Attn:  S.  0.  Dean  (3  copies) 

G.  Hess 

R.  L.  Hlrsth  (3  copies) 
B.  Miller  (3  copies) 

R.  Blacken 

Lawrence  Livermore  Laboratory 
University  of  California 
Livermore,  California  94551 

Attn:  T.  K.  Fowler 
D.  Baldwin 

H.  Berk 

J.  Killeen 
B.  McNamara 


J.  Rome 
G.  Bateman 
C.  Beasley 
R.  Dory 
C.  Cruoe 


R.  Price  (3  copies) 
V.  Sadovskl 
D.  Pries ter 
A.  Sleeper 
A.  Davies 


R.  Post 
J.  Byers 
W.  Kruer 
E.  Valeo 
B.  Langdon 


Massachusetts  Institute  of  Technology 
77  Massachusetts  Avenue 
Cambridge,  Mass.  02139 

ATTN:  B.  Coppi  (Riysics)  D, 

T.  H.  Dupree  (Njc.  Eng)  J. 

A.  Bers  (E.  E. ) R. 


U.  C.  L.  A. 

Los  Angeles,  California  90024 


Attn:  J.  Dawson 
B.  Fried 
Y.  C.  Lee 
F.  F.  Chen  (E.  E.) 


D.  Slgmar  (Nuc.  Eng.) 
J.  McCune  (A.  E.) 

R.  Parker  (E.  E. ) 


1 


/ 


Princeton  University 

Plasma  Physics  Laboratory 

Princeton,  New  Jels«y  08540 

Attn:  M.  B.  Gottlieb 
R.  Ellis 
Z.  Friesian 
H.  Furth 
D.  Grove 
D.  Haade 
JL  Mills 
P.  Rutherford 
T.  Stiac 

Los  Alamos  Scientific  Laboratory 

P.  0.  Box  1663 

Los  Alamos,  New  Mexico  87344 

Attn:  P.  L.  Kibe 
D.  Forslund 
C.  Nielson 
J.  Friedberg 
R.  Morse 
J.  Kindell 


D.  Jassby 
V.  Tang 
F.  Perkins 
J.  Johnson 
X.  Kulsrud 

B.  Okuda 
T.  K.  Chu 

C.  Obernan 
J.  M.  Greene 


B.  Godfrey 
R.  Mason 
W.  Ellis 
E.  Llndman 
R.  R.  Lewis 
J.  U.  Brackbill 


University  of  Texas 
Austin,  Texas  78712 


Attn:  D.  Ross 

W.  Drummond 
A.  Ware 
L.  Sloan 


F.  Hinton 
W.  Horton 
R.  Hazel Cine 


9.  University  of  California 
Berkeley,  California  94720 

Attn:  Prof.  A.  Kaufman 

Prof.  C.  K.  Blrdsall 


10.  Bell  Telephone  Laboratories 
Shippany,  New  Jersey  07981 

Attn:  S.  Buchsbaum 
A.  Hasegawa 


..J.  ■ ....  ■■  J-  T-aiRy  . .1 


* 


2 


11.  National  Science  Foundation 
Washington,  D.  C.  20550 


Attn:  E«  C.  Creutz 
R.  Sinclair 

12.  General  Atonic  Co. 

P.  0.  Box  81608 

San  Diego,  California  92138 

Attn:  T.  Ohkava  M.  S.  Chu 

D.  Dobrott 
G.  Guest 

13.  Weselnghouse  Electric  Conpany 

Breeder  Reactor  Division 

P.  0.  Box  355 

Pittsburgh,  Pennsylvania  15230 

Attn:  Z.  Shapiro 
D.  Klein 
R.  Rose 
T.  Varljen 

14.  University  of  Wisconsin 

Madison,  Wisconsin  53706 

Attn:  R.  Conn 

G.  Kulclnski 


15.  New  York  University 

Ccurant  Institute  of  Math.  Sci. 
251  Mercer  St. 

New  York,  N.  Y.  10012 


Attn:  H.  Grad 

W.  Crossoar. 

H.  Weitzner 
J.  Tataronis 

16.  Harwell  Laboratories,  Inc. 
9244  Balboa  Ave. 

San  Diego,  CA  92123 


Attn:  A.  Kolb 
V.  Fargo 
N.  Rostoker 


■ vt„ 


m 


> 

A 


I i 
I ' 

m.  .*‘j 

I 


i 


17.  Institute  for  Advanced  Study 
Olden  Lane 

Princeton,  New  Jersey  03340 
Attn:  M.  Rosenbluth 

18.  Ba telle 

Pacific  Northwest  Lab. 

P«  0.  Box  999 

Richland,  Washington  99352 

Attn:  L.  Schmidt 

19.  Foraen  Incorporated 
777-106  Avenue,  H.  W. 

Bellevue,  Washington  98004 

Attn:  H.  For sen 

20.  University  of  Maryland 
College  Park,  Maryland  20742 

Attn:  P.  C.  Llewer  C.  S.  Liu  D.  Tidmaa 

Prof.  R.  Davidson  H.  Grlea 

21.  Stevens  Institute  of  Technology 
Hoboken,  Hew  Jersey  07030 

Attn:  George  Schmidt 
B.  Rosen 

22.  College  of  William  and  llary 
Williamsburg,  Virginia  23185 

Attn:  S.  P.  Gary 

Fred  Crown field 

23.  Culham  Laboratory 
Abingdon,  Berkshire,  England 

Attn:  Chris  Lasrcore  Davies 
Ian  Cook 
Keith  Roberts 
Jes  Christiansen 
J.  B.  Taylor 
B.  Peas 


M.  Hughes 
J.  A.  Wesson 
T.  Stringer 
J.  W.  M.  Paul 
A.  Gibson 


24.  Sand la  Laboratory 
Albuquerque,  New  Mexico  87115 

Aetn:  Thomas  Wright 
J.  Freeman 

25.  University  of  Rochester 
Rochester,  New  York  14627 

Attn:  Moishe  Lubln 

E.  B.  Goldman 

F.  J.  Cat to 

26‘  ^Pi‘Qck  Iaacit'lte  i^r  Plasma  Physics 
8046  Carching  bel  Munchen 
West  Germany 


Attn:  D.  Biskamp  y 

Norst  Packer  C 

W.  Schneider  D*  P 

27.  Physics  International  Company 
San  Leandro,  California  94577 

Attn:  V.  Bailey 

28.  Science  Applications,  Inc. 

Laboratory  for  Applied  Plasma  Physics 
Research  Staff 

La  Jolla,  California  92037 


K.  Von  Hag enow 
Chodura 
D.  Pfirsch  ‘ 


Attn:  H.  Krall  c,  Wagner 

R.  Shanny  H.  Klein 

J.  McBride  N.  Byrne 

29.  University  of  Tokyo 
Bunkyo-ku 
Tokyo,  Japan  113 

Attn:  Solchl  Yoshikawa 

30.  Japan  Atomic  Energy  Research  Institute 
Tokai-!!ura 
Ibarakl-Ken,  Japan 

Attn:  N.  Fujisawa 


31.  Association  Euratoa-CEA  Serv  La  Fusion 
Dept,  of  Physics  du  Plasna  at  de  la 

Fusion  Controlaa 
Center  d’ Etudes  Kuclealres 
Bolte  Postale  Ko  6 
92260  Fonteney-Aux-Rose 
France 

Attn:  Dr.  Tachon 
Dr.  Kercler 
R.  B.  Paris 

32.  Australian  Rational  University 
Canberra  A.  C.  T.  Australia 


Attn:  J.  D.  Strachan 

R.  L.  Dewar 

33.  Kurchatov  Institute 
Moscow,  U.  S.  S.  R. 

Attn:  B.  B.  Kadomtsev 
E.  P.  Velikov 

34.  Cornell  University 
Ithaca,  Hew  York  14850 

Attn:  Prof.  R.  N.  Sudan  C&PP*  Phys.) 

Prof.  E.  Ott  (E.  E.) 

S.  V.  Lovelace 

35.  Kaval  Research  Laboratory 
Washington,  D.  C.  20375 

Attn:  Code  7700  (Tens  Mason)  - 25  copies 

Code  7750  - 75  copies 


I 'I 


AAVN  MU  tO 


PENALTY  FOR  PRIVATE  USE.  800 


